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Abstract 

We develop a scheme to distill entanglement from bipartite Fermionic 
systems in an arbitrary quasifree state. It can be applied if either one 
system containing infinite one-copy entanglement is available or if an ar- 
bitrary amount of equally prepared systems can be used. We show that 
the efficiency of the proposed scheme is in general very good and in some 
cases even optimal. Furthermore we apply it to Fermions hopping on an 
infinite lattice and demonstrate in this context that an efficient numerical 
analysis is possible for more than 10 6 lattice sites. 



1 Introduction 

An important class of Fermionic states are the Gaussian ones. They describe 
equilibrium states of quasifree spin chains, such as non- interacting electrons |16j . 
Since they are completely characterized by their two-point correlations, they are 
under good control also in case of large particle number. For an example of their 
standard description based on a fixed basis in the Fermionic Fock space, see [7]. 

Defining entanglement is not straightforward in Fermionic systems as 
Fermions are indistinguishable and described by antisymmetric tensor prod- 
ucts. There are several ways around this problem. E.g. in |llj antisymmetrized 
states with definite particle number are studied while in [T| the focus is on the 
lattice spin system admitting Fermionic description. The conceptionally clear- 
est approach is, however, to base the description of subsystems not on tensor 
products of Hilbert spaces, but on the notion of local observables [5]. This ap- 
proach is successfully applied to the study of entanglement of systems with 
infinite degrees of freedom |20[ [P3] and for the analysis of separable [T7] and 
maximally entangled [TB] states of Fermionic systems. Furthermore, the conser- 
vation of local parity superselection rule allows for different possible definitions 
for separability [5J (consult that article also for more literature on Fermionic 
entanglement). 

In this paper we will take an operational point of view. Instead of asking 
how much entanglement is contained in a given bipartite, Fermionic system, we 
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will look for explicit distillation protocols, i.e. procedures to generate (in terms 
of LOCC), from a number of Fermions, pairs of distinguishable particles in a 
maximally entangled state. The asymptotic rate of such a protocol can then 
be regarded as a measure for the entanglement contained in the original state. 
The advantage of this approach is that the only concept which needs a (slight) 
generalization is that of local operations (of bipartite Fermionic systems and 
from a Fermionic to an ordinary bipartite system). The latter, however, can be 
easily based on the idea of local observables mentioned in the last paragraph; 
cf. |12) for a more complete discussion. 

Following this idea, we present a family of protocols which is particularly 
adopted towards distillation from quasifree (i.e. Gaussian) states. The general 
structure involves a two step procedure: First we generate (with a certain success 
probability p) bipartite d-level systems (where the local dimension can be chosen 
within a certain range) in an isotropic state a. In particular for very large 
systems (cf. the discussion of free Fermions on an infinite lattice in Section [5]) 
the entanglement fidelity of a can be already very close to one such that no more 
steps are necessary. If this is not the case but a is distillable we can continue 
with standard techniques like the recurrence method or hashing (cf. [5T] and 
the references therein). 

Each protocol in this family admits a very easy parametrization in terms of 
two operators (a projection D and a partial isometry V; cf. Section!?]) on the one- 
particle reference space K, rather than the corresponding Fermionic Fockspace 
W. This implies in particular that we can express the success probability p and 
the fidelity / of the isotropic output states a (from which the overall rate of 
the protocol can be calculated with known formulas, once we have decided how 
to process the resulting ci-level systems) in terms of D and V. The dimension 
of JC (i.e. the space on which D and V operate) grows, however, only linearly 
in the system size (i.e. the number of independent modes) while dim'H grows 
exponentially. Hence we get a very efficient way to discuss the entanglement 
content of quasifree Fermionic systems which is applicable to, literally, millions 
of particles. This is explicitly demonstrated in Sectional 

The given class of protocols is still large and for a given quasifree state ps 
(which is completely described by its two-point correlation matrix S) most of 
them lead to very poor results. Therefore, in Section [5] we present a scheme to 
derive the operators D and V from S. It is based on the observation that to any 
generic, quasifree state ps we can associate a quasifree, maximally entangled 
state in a natural way (this resembles a little bit the Schmidt decomposition, 
but note that pureness of ps is not required). For a large subset of quasifree 
states (whose dimension grows quadratically in the system size) this particular 
protocol optimizes the product of success probability p and fidelity / (we will 
discuss in Section [5] why this is a good figure of merit and, actually, related to 
the overall distillation rate) . 

The presented scheme is discussed in terms of two classes of examples: Small 
systems consisting of two or four modes and free Fermions, hopping on a one- 
dimensional lattice. The numerical efficiency mentioned above allows us in the 
latter case to treat large particle numbers and this leads to new insights about 
the entanglement properties of this particular system, which are interesting in 
their own right. 

The paper is organized as follows: In Section [5] we will start with some 
mathematical preliminaries about Fermionic systems. The approach we present 
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here is based on the "selfdual formalism" [31 [5] and might be unfamiliar for most 
readers (although it is not particularly new). For the purposes of this paper 
it is, however, the most natural formalism. In Section [3] we apply it to discuss 
some aspects of bipartite Fermionic systems and in Section 0] we present the 
family of protocols mentioned above. The question of optimality is then treated 
in Section [SJ The rest of the paper is finally devoted to examples: two and four 
modes in Section [5] and Id lattice systems in Sections [7] and [SJ Several technical 
aspects (including, in particular, proofs) are postponed to the appendix. 

2 Mathematical preliminaries 

To introduce some terminology and notations, let us start with some technical 
remarks about Fermionic systems. All the material we are going to present here 
can be found in the literature, [2] are some standard references. 

A Fermionic system consisting of n modes can be described by smeared 
Majorana fields 

K, {n) 3x^B{x) eB{H (n) ). (1) 

The system Hilbert space ?{W is 2™ dimensional and can be realized as the 
Fermionic Fockspace over C™. Therefore we will refer to it occasionally as the 
Fermionic Fock space of n modes although the precise form of W'"' is (apart 
from its dimension) not really important. The Hilbert space - in the follow- 
ing called reference space - is 2n dimensional and equipped with an antilinear 
involution (called complex conjugation) J : /C( n ) — > KP 1 '. In other words J 
satisfies 

J(x + Ay) = J{x) + XJ(y), J 2 = 1 (2) 

for all x, y e K,^ and all A e C. The typical choice for and J is C 2 " 

equipped with the ordinary conjugation in the canonical basis. As with H^ n \ 
however, the explicit realization of [K^ n ',J) is not important. K.^ contains a 
distinguished real subspace given by 

/C (n) \Jx = x}c JC {n) . (3) 

Lots of structures we will encounter in the following are actually associated to 
JQg rather than to lC^ n \ We will call in particular an orthonormal basis ej, 
j = 1, . . . , 2n of /Cg™' (which is of course an orthonormal basis of KS n \ too) a 
real basis. 

The operators B{x) are (complex) linear in x € JQ- n ' satisfy the canonical 
anticommutation relations (CAR) in the form 

{B(x),B(y)} = (Jx,y)l, B(x)* = B(Jx) (4) 

and they act irreducibly on H^ n \ i.e. 

[A, B{x)] = Vx e JC {n) =*> A = Al. (5) 

These conditions fix the B{x) up to unitary equivalence, which is the reason 
why we are not interested in their explicit form. They can be constructed easily 
in terms of ordinary creation and annihilation operators. The details are shown 
in Appendix [Al 
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Let us consider now states of the system. To any density operator p on T-L^ 
we can associate a covariance operator S G B(K.^ n ') by 

tv(pB(x)B(y)) = (Jx,Sy) Vx,y € K. (6) 

S is selfadjoint and satisfies 

< S < 1 and JSJ = 1-5. (7) 

A state ps is called quasifree if it is uniquely characterized by S and the condi- 
tions 

ti(p s B(xi)---B(x 2l+1 ))=0 (8) 

l 

tr(p s B(xx) ■ ■ ■ B(x 2 i)) = ^sign(p) Y[( Jx P (2j-i), Sx P (2j)), (9) 

i=i 

which have to hold for all I € N and Xk & IC^ . The sum in © is taken over all 
permutations p satisfying 

p(l) <p(3) < ... <p(2Z-l), p(2j - 1) < p(2j) (10) 

and sign(p) is the signature of p. 

A pure state ipE G is quasifree (and then called a focfc state) iff its 

covariance operator is a projection. According to Equation it has to satisfy 
JEJ = 1 — E. Each projection with this property is called a basis projection. 

Consider now a unitary operator R on. K. satisfying [J, R] =0. It leaves the 
real subspace /Cr invariant; i.e. it is a rea? orthogonal transformation. R gives 
rise to a new set of operators by Br(x) = B(Rx). It is easy to see that Br 
is complex linear and satisfies (jU) and ([S]). Hence the fields B(x) and Br(x) 
are unitarily equivalent and describe effectively the same physical system. More 
precisely there is a unitary T(R) on HS n \ called the Bogolubov transformation 
of R satisfying 

a R {B{xj) = T(R)B(x)T(R)* = B{Rx). (11) 

The Bogolubov automorphism cir of B{Jv- n ') is uniquely determined by this 
condition, the Bogolubov transformation is only fixed up to a phase. 

The most important special case arises with R = — 1. The corresponding 
automorphism O = a_j is called the parity automorphism. According to pip it 
is characterized by 

Q(B(x)) = 6B{x)6 = -B(x). (12) 

The associated Bogolubov transformation r(— 1) = 9 can be chosen selfadjoint 
and is then called parity operator. It is fixed by (|12[) up to a sign and given in 
terms of a real basis e a , a = 1, . . . , In of KS n ' by (cf. Appendix IB"]) 

6 = 2 n i n B( ei )---B(e 2n ). (13) 

If we change the basis in terms of an orientation preserving, real, orthogonal 
transformation, 9 remains invariant. If we change the orientation, the operator 
9 changes the sign. In the following we will assume that a particular orientation 
is chosen. 
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The parity automorphism gives rise to the distinction of even and odd op- 
erators: A G B{U (n) ) is called even if 0(A) = A and odd if Q(A) = -A. Only 
selfadjoint, even operators can be regarded as observables. Similarly, a com- 
pletely positive (cp) map T : B{W^) — > B(H.( n ^) is an operation only if it 
commutes with 0, and only its action on even elements is relevant. In other 
words, if a second cp map T\ satisfies T(A) = T\{A) for all A with <d(A) = A it 
describes the same operation. 

Using the spectral decomposition 9 = P + — P_ of the parity operator we can 
decompose the Hilbert space VS n ^ into an even and an odd part: 

n = n { + } ffi n { l l) , h { ± i) = p±n {n) , e = p+ - p_ . (w) 

An operator A = B(H^) is even iff it is of the form A = A + ® A_ with 
A± eB(n^ y ). 

A Fermionic subsystem of our given system consists of a number of modes 
which are distinguished by a certain physical property like position in space. 
Mathematically it is described by a projection D satisfying 

DeB(IC), D 2 = D, D*=D [D,J]=0. (15) 

The last condition implies that D can be restricted to the real subspace K,^ of 
Therefore we will call it a real projection. D projects onto the subspace DK, 
of JC( n ' containing the modes belonging to the subsystem. The corresponding 
Majorana operators does not act any longer irreducibly on "H^. Instead, we 
can find a unitary 

U D :H {n) -^n {l) ®% {n ' l \ l=tr(D)<n (16) 

which satisfies and is (up to a phase) uniquely characterized by 

U BMW \ B{X) ® 1 X G DK(n) ~ ^ h7) 
UdB{x)Ud ~ \9 ® B(x) *e(I-i?)JC(»)sJC(»-'). (17) 

Here n (l) and H (n ~ l) denotes the Fockspace of I = tr(D) and n - I = tr(l - D) 
Fermionic modes. Similarly K,^ and JC( n ~ 1 ^ are the corresponding reference 
spaces. 

Dropping the complementary subsystem is described by the operation (in 
the Heisenberg picture) 

B{H {1) ) 3i4 Ajj(A) = U* D A ® 1U D G B{U {n) ). (18) 

It is easy to see that even operators on are mapped to even operators on 
n (n) . Hence the map describes an operation in the sense discussed above. 



3 Bipartite systems and entanglement 

Let us consider now a bipartite systems consisting of 2m modes shared by 
Alice and Bob. Hence n — 2m and there is a distinguished n-dimensional, real 
projection I a G B(K,( 2m ' > ) which defines the Alice subsystems. Likewise 1b = 
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I — A defines the Bob subsystem, and the Hilbert space K.^ 2 " 1 ^ decomposes into 
a direct sum 

£(3m) = K W (BK W > J = J A(S J Bj = I A/B K,^ m \ J A/B = I A/B JI A/B . 

(19) 

According to Equation (fT7|) we can identify the Hilbert space "H^ 2 ™) via the 
unitary U A — Uj A with W m ' W m ' such that ordinary entanglement theory 
applies. There are, however, two important caveats. 

According to (jTTJ) a Majorana operator B(x) belonging to Bob's subsystem, 
i.e. x G I B IC {2m) is of the form U A B{x)U* A =6® B(x). Hence Alice's Hilbert 
space is not invariant under the action of Bob's operator^. This problem can be 
rectified if we take the remark seriously that only even operators are physically 
relevant. Since an operator Q £ B^MP 1 "^) is even iff it can be written as an even 
polynomial in the B(x), it follows immediately that an even operator belonging 
to the Alice (Bob) subsystem operates trivially on the Bob (Alice) Hilbert space. 
Or in other words the observable algebras 

A = span{S(a; 1 ) • ■ ■ B(x 2 i) \ x u x 2 i G I A K, I G N} (20) 
B = span{B(a; 1 ) ■ • • B(x 2l ) | x u x 2 i G I B K, I G N} (21) 

associated to Alice and Bob respectively are of the form 

U A AU* A = (S(?4 ro) ) 8 6(H W )) ® I c B{U {m) ) ® I (22) 

U A BU* A = 1® (B{U { + l) )@B{U {m) )) c l®B{H (m} ). (23) 

But now a second problem arises, since both parties admit observables - the 
local parities 9 (g> 1 and 1 ® — which can be measured without disturbing the 
system. Therefore we have to deal with entanglement theory at the presence of 
local superselection rules. A detailed discussion of this subject can be found in 
|19j . For us only a few special topics are relevant. 

Let us start with a quasifree state ps with covariance operator S. The sub- 
system projections I A and I B give rise to the operators Sjk = IjSIk, j, k — A, B. 
The properties of S imply that 

X s = -i(2S A A-t), Zs = -i(2S B B-t), Y s = -2iS AB = 2iS* BA (24) 

are real operators and that X$ and Z$ are antisymmetric. 

Many entanglement properties of ps can be stated in terms of Xs, Y$ and Z$. 
Important for us are maximally entangled, quasifree states which are according 
to [TB] characterized by the condition 

X s = Z s = 0, Y S Y*=I A , Y*Y S = I B . (25) 

In other words Ys is a partial isometry with I A as source and I B as its target 
projection. If we introduce a real basis e a , a = 1, ...,4m of K,^ 2 ™^ which is 

x We will do this in the following without further reference, as long as confusion can be 
avoided. 

2 Semingly, a similar problem does not arise with Alice's operators. This is, however, only 
an artifact of the special representation given by the unitary U A ■ We could equally well choose 
Ub = Ui B to exchange the roles of Alice and Bob. The crucial point is that Majorana operators 

B(x A / B ) with x A j B £ ^C^T/g do not act independently. 
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adopted to the Alice/Bob split, i.e. 

ei,...,e 2m eI A rt 2m \ e 2m+1 ,...,e 4m e/ B /C( 2m ) (26) 
we get the matrix 

Y jk = (ej,Y s e 2m+ k), j,k = l,...,2m (27) 

which is orthogonal. This gives us a parametrization of the set of all quasifree, 
maximally entangled states in terms of the group 0(2m), and it implies that 
there are two possible orientations. This is connected to the global sign of 9 
which can be fixed by the condition tr(9ps) > (provided this expectation 
value is not zero; note in this context that ti(pgO) — ±1 only, occurs iff ps is 
pure, i.e. a Fock state). 

According to the definition each quasifree state is an even state, i.e. 
tr(psA) = for each odd operator A on 'H^ 2m \ This implies P+psP- = 
P-PsP+ — 0, and P+i/je = or P-ipE = if p = \iPe)(iPe\ is pure. Note 
that the roles of P + and P~ can be exchanged by switching the orientation 
(and therefore the sign of 9; cf. the discussion of Equation ([T5)l ). If we fix the 
sign of 6, as stated in the last paragraph, by the condition that (tpEtOtpE) = 1 
holds we get P~ipE — 0, which should hold in the following. 

Now consider the decomposition #( 2m ) = (g)#( m ) of the global parity into 
a product of local ones. It implies a likewise decomposition of the projections 
P± as 

P+ = P++ + P— , P- = P+- +P-+, Pjk = Pj ®Pk, j,k = ± (28) 
and therefore P-ipE = implies 
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4> E = 4= + fa,-) , ^em g n { k m) ® H[ m} , k = ± (29) 



and ih_E is maximally entangled if ip E ,± ^ re maximally entangled in the usual 
sensqj. 

An important point is now that a relative phase between ipE,+ an d ipE,— can 
not be determined by local measurements by Alice and Bob. In other words ipE 
and 

i> = -±= {i> E ,+ + e- ia i> B ,-) (30) 

are completely equivalent, as long as only LOCC operations and measurement 
are allowed. The only quasifree choice in this family arises for a = n. The 
corresponding basis projection E differs from E only by the sign of the off- 
diagonal block Eab> i- e - 

Saa = Saa, Sbb = Sbb, Sab = —Sab- (31) 

If ipE is maximally entangled is maximally entangled as well and it represents 
the same orientation, i.e. {ips,,6tp^,) — {ipE,9ipE)- 

Finally, let us have a short look at LOCC. As with entanglement, the usual 
concepts apply to Fermions if we use the tensor product decomposition of %( 2m ) 



3 Note, however, that not all maximally entangled states can arise here, since ipE is quasifree 
by assumption. 
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given by Ua and take into account that only operators in the tensor product 
A® B are physically relevant. We are not giving a full discussion here, but refer 
the reader to [T2|. Instead we will only have a short look on those operations 
which are relevant for our purposes. 

Our first example is "dropping a subsystem" as described in Equation (|18[). 
If the subsystem projection D commutes with I a we can decompose D as D = 
Da © Db with Da — IaD and Db = IbD. The overall operation can 
therefore be written as Ax> B . The final system consists of tr(Z? J 4)+tr(Ds) 

Fermionic modes. 

Now consider the projection operators Pjk, j, k = ± introduced in Equation 
(|28[), They define a von Neumann-Liiders instrument which we will call in the 

following a joint parity measurement. The possible values are H — h, H — , — h, , 

and if the system is before the measurement in the state p we get 

p> k = tv(P^p) and p> k = ^— (32) 

pj 

for the probability p jfe to get the outcome jk and for the corresponding posteriori 
state pP k . Each of the subchannels p H> PjkpPjk is obviously a local operation. 
They produce bipartite systems, which are described by the Hilbert spaces 

?4 m) <8> n{ m) = C d <g> C d , d = 2'"- 1 . (33) 

Hence we get a pair of (distinguishable) d— level systems. For later reference let 
us introduce as well the probability 

p = P+++P— = tr(p(P ++ + P__)) = tr( P P+) (34) 

to get the same parity on both sides. Since P++ + P__ = P+ this definition 
depends on the choice of an orientation (i.e. the sign of 9). 

An LOCC operation which produces again a Fermionic system (but in con- 
trast to Ac of the same type) is twirling. We choose a quasifree, maximally en- 
tangled state pe = \iPe)(iPe\ and average over the group G of all local unitaries 
U = U A ® U B on -H (2m) = U {m) ® H (m) such that 1^-4^ = .4, f7 s BL/£ = 6, 
and f/V'E = "0b holds. This leads to a state 

T E p= ( UpU*dU (35) 
Jg 

= \+\i> E }(<lpE\ + \-\Tp E ){lp E \ + M+(P ++ + P—)+ P-(P + - + P- + ) 

(36) 



with 



p A+ + A_ , 2 l-p _ 



■ M+rf > = M-rf , (37) 

(ip E ,atp E ) = A++ /x+, (ipz,o-i/jz) = A_ (38) 
where p is the probability from (|3"4"|) and <i is given as above by: 2 m_1 . 



4 A family of distillation protocols 

Let us consider now a large number of bipartite Fermionic systems, each consist- 
ing of 2L mode^f], which are macroscopically distinguishable (e.g. metallic wires, 

4 Basically, we could also use an infinite number of modes at this point, but some of the 
statements made in the last two sections are not valid in this case. 
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each containing an electron gas), and each prepared in the quasifre^l state ps 
(e.g. a thermal equilibrium state or a ground state of the electron gas). To distill 
from these systems maximally entangled pairs of d-level systems for some d £ N 
we can proceed as follows (cf. |12|): 

• Step 1. Choose an integer 1 < m < L and a real projection D £ B{JC^) 
which commutes with the projection 1a onto Alice's modes, and sat- 
isfies \x{T>a) — tr(Ds) — m. It defines a bipartite subsystem of our 
original system and we can drop the remaining modes by applying the 
channel Ac denned in Equation (|T51) . We get a Fermionic system con- 
sisting of 2m modes in the quasifree state pdsd with covariance matrix 
DSD £ S(/C (m) ) B{DJC^). The purpose of this step is to discard low or 
unentangled parts of our original system. The correct choice of D requires 
knowledge about the state p$. 

• Step 2. Now choose a partial isometry V £ B(K.( 2m ^), satisfying VV* = 
Da and V*V = Db- It defines a maximally entangled, quasifree state 
ipE £ T-L^ 2m ^ of the reduced system. It has to satisfy 

P 

{^e,Pdsd^e) + {iPe,P^e) > ( 39 ) 

otherwise the protocol would fail. Now apply the corresponding twirling 
operation (|3"Tj) to get the state Te(pdsd) from Equation 

• Step 3. Make a joint parity measurement. This leads to a pair of distin- 
guishable systems described by the tensor product of two d-dimensional 
Hilbert spaces. 

If the outcome is H — h or (which happens with probability p given in 

Equation (|34p) this leads to the isotropic state er ++ or a given by 



^ = ^-^IfeKfel + -fei^, (40) 

zp p 

/( m ) <i/( m ) 



where the maximally entangled states ipE,k £ H\. ® n k are defined by 
the decomposition (cf. Equation (j2"9"]) ) 

l/j E = i(^E,+ +^ B -) (41) 

The fidelity / = (ipE,±, & ±± ' 1 Pe,±) of a ±± is given by 

j = (^e,Ps^e) + (iPeiPs^e) ^ 
P 

such that Equation ([55)) becomes / > h. Hence is distillable. 

Technically, a ++ and a are denned on different Hilbert spaces. Both, 
however, differ only by the value of the local parity operators 9 <S> TL and 
I <£) 8. If we ignore the value of the parities (and only look whether they 



5 Basically the protocols we are going to describe work for non-quasifree states too, but 
they are probably not very good. Furthermore, all non-trivial statements we will make in this 
paper are related to quasifree states. 



9 



coincide or not) we can identify ® U% n) and U { l n) ® H M such that 
4>e,+ and ipE,— , and tr ++ and a become identical. 

If the outcome of the measurement is H — or — + we get a totally mixed 
state which has to be discarded. 

• Step 4. Take the next system and restart at Step 1 (using the same D 
and the same V) . Repeat this procedure until enough systems in the state 
<7 ±=t are available for standard distillation techniques like hashing [3T] and 
proceed accordingly. 

Hence, for each integer n describing the size of the input systems, we have 
defined a family of distillation protocols which is parametrized by the integer 
m, the real projection D and the partial isometry V. A useful choice of all three 
parameters requires knowledge about the covariance matrix S of the input state 
ps- The quality of the outputs is then measured by the probability p(S, m, D) 
from Equation (|54"1) and the fidelity f(m, D, V, S) given in (|4"2"|) . Both quantities 
can be calculated easily according to the following two theorems. 

Theorem 4.1 For each quasifree state ps of 2m Fermionic modes with co- 
variance matrix S the probability p — tv(psP+) to get the result + in a parity 
measurement ( cf. Equation \3J$ ) is given by 

l + (-4rPf(-i(S-I/2)) 
P= ^ — i (43) 

where Pf denotes the PfafharE 

Proof. See Appendix [B] □ 

Theorem 4.2 The fidelity of a quasifree state ps with covariance matrix S with 
respect to a quasifree, pure state tpE with covariance matrix E is given by 

(il>E,Ps1>E)=~Pf(-i(Z-S-E)). (44) 

Proof. See [IS]. □ 

Hence the two quantities we are interested in become 

P (S,m,D)= 1 + ^ mFf i DSD - 1 ^ (45) 



tra n,/i 1 wf d aXsD a D a Y s D b + V \ 
f(S, m, D,V) = ¥ ^Pf{ _ DbY * Da _ y DbZsDb ) + 

J_ / D A X S D B D A YD B -V\ 

22m rl ^ -D B Y*D A + V D B Z S D B ) ' yw> 

where we have used the fact that Xdsd — P>aX 8 Da and similar relations for 
Ys and Zs hold. 

6 Note that the Pfaffian can be defined in a basis free way for any antisymmetric operator on 
a real, even-dimensional, oriented Hilbert space. The —i factor under the Pfaffian in Equation 
1 143 I I is necessary to get a real (rather than purely imaginary) operator. 
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Before we consider the question how to choose the projection D and the 
partial isometry V , let us have a short look on a variant of the protocol given 
above. We can skip the twirl in Step 2 and perform the joint parity measurement 
on the system in the state pdsd- If the outcome is H — or — h we drop the 
systems as before, and in the other case we get with probability p systems in 
the state 

P= —Pdsd + —Pdsd-. ( 47 ) 

where we have identified (as above) the Hilbert spaces ®H+ and ® 
H_ such that ipE,+ and ipE,- coincide. If we twirl now with the symmetry 
group of ipE,± we get the same output states a as in Step 3 above. In other 
words, nothing has changed, apart from a reordering of steps. We will come 
back to this point, however, in Section [5] where we will compare our protocol to 
a slightly different one. 



5 The optimal protocol 

To distill entanglement from a given quasifree state ps with a protocol from the 
family just introduced, the projection D and the partial isometry V have to be 
chosen carefully, otherwise we do not get anything. The purpose of this section 
is to show that there is a canonical choice, which provides good or even optimal 
results. 

The crucial point is that any S provides, by the polar decomposition 

Y s = Vs\Y s \ (48) 

of the off-diagonal block Y — — 2iSAB (cf. Equation (|24])1. a partial isometry 
which is uniquely determined by the condition 

Ran^s = Ran|y s | C Ranis and kerVs = ker|r 5 | C RanI A - (49) 

Hence, any covariance matrix S satisfying Ran|Ys| = Ib and ker|Ys| = I a 
(which holds iff \Ys\ has maximal rank) defines a unique quasifree, maximally 
entangled state via the partial isometry Vs- This observation immediately sug- 
gests the following procedure to select D and V: 

• Subsystem projection D. Consider an eigenbasis e*, k = 1, ...,4L of 

\Yg\ such that the corresponding eigenvalues are arranged in decreasing 
order, i.e. 

\Y s \e k = X k e k , Ai > A 2 > ••• > \ iL = 0. (50) 

Note that the kernel of Ys is by construction at least 2L dimensional. 
Hence X2L+1 = ■ ■ ■ = X^l = 0. Now choose 2 < m < n small enough, but 
at least such that Afc > holds for all k < 2m (we rediscuss strategies for 
the choice of m below). With this preparation we can define D = Da®£>b 

by 

2m 

D s = ^|e fe )(e fc |, Da^VsDbVI (51) 
fc=i 

This construction is in most cases basis-independent {Db is just a spec- 
tral projection of |Ys|). The only exception arises if the eigenvalue X m 
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is degenerate. In this case L>b depends on the choice of a basis in the 
corresponding eingenspace. 

• Partial isometry V. Now consider the restricted state defined by the 
covariance matrix DSD, and choose V such that 

Y bsb = v\ Y bsb\ with Y bsb = -2i{bsb) AB = b A Y s b B , (52) 

i.e. the maximally entangled state defined by V is the "natural" one con- 
sidered above. 

We do not know yet, whether this choice really leads to a good distillation 
rate, or whether it is even optimal. If we choose m — 2 in Step 1 (i.e. the out- 
put systems are qubit pairs) and hashing as the standard distillation technique 
needed in Step 4 we get (where S(a ±± ) denotes the entropy of the isotropic 
state (7±± from Equation ([1D])1 

R = pmax{0, (l - 5(a ±:t ))} (53) 
= pmax{0, (1 + /log 2 (/) + (1 - /) log 2 (l - /) - (1 - /) log 2 (3))}. (54) 

for the overall distillation rate (cf |21) for the rate of hashing in any prime 
dimension). There are two problems with this quantity. First of all if S(a ) > 
1, hashing does not lead to a nonvanishing rate, although other strategies (like 
applying the recurrence method first) might be more successful. Optimizing R 
is pointless in these cases. Second, even if R > holds it is a very difficult 
function of D and V (cf. the expressions for / and apply Theorems 14.11 and 14 . 2[) 
and optimization is therefore a hard task. 

From a more general point of view, however, R is just a product of p and a 
quantity which is monotone in /. This property it shares with 

Pf = {tpEiPs^E) + {^e,PSiPb), (55) 

which is independent of a special distillation strategy in Step 4 and easier to 
optimize. As the rate R it avoids strategies where the fidelity / is increased at 
the cost of the probability p (which is possible, according to Equation (02)). It 
shares with other fidelity quantities the drawback that it can captures effects 
which are not related to entanglement at all (e.g. local manipulation). In our 
case, however, and in particular for large /, it provides a reasonably good figure 
of merit, which allows - at least for a large subclass of quasifree states ps - the 
quantities D and V as optimal solutions. More precisely we have the following 
theorem (the proof is given in the appendix) . 

Theorem 5.1 Consider 2L Fermionic modes in a quasifree state ps, an integer 
2 < m < L, operators D, V , and the quantities p(m, D, S) and /(to, D, V, S) 
from Equations J^<5| ) and f^6] j. If Xs = or Zg = (cf. the notation from 
Equation \2J$ ), the product pf satisfies 

2m 1 4. \ 2m i _ \ 
p(m, D, S)f(m, D, V, S) < p(m, D, S)f(m, D, V, S) = JJ + [] —f 1 , 

k=l k=l 

(56) 

where b, V are defined above, and the Afe are the 2m highest eigenvalues of \Yg\ 
(cf. Equation 15(A I). 



12 



Proof. See Appendix [Cl 



□ 



Up to now we have kept the number m (i.e. the modes left after Step 1) 
fixed. The reason is that a good choice for m depends much more than the 
choice of D and V on the initial state ps and the distillation strategy used in 
Step 4. The crucial question is up to which degree it can be more beneficial to 
increase the dimension d of the output system at the expense of some fidelity (cf. 
the corresponding discussion in |21|). If, however, we are only going for fidelity 
m = 2 is the best choice. 

Corollary 5.2 Consider the same assumptions as in Theorem \ 5.1l Then we 
have 

p(2, D, S)f{2, D, V, S) > p(m, D, S)f(m, D, V, S). (57) 

Proof. See Appendix [Cj □ 

The condition X$ = or Z$ — used so far is restrictive, but the corre- 
sponding class of states is still fairly big and grows quadratically in L. Further- 
more it seems to affect only local correlations (i.e. only between Alice's modes or 
only between Bob's modes). Therefore we conjecture that the protocol given by 
the operators D, V presented above is generically a good if not the best choice 
within the family introduced in Section 2) We will refer to it in the following 
as the optimal protocol, even if the state to distill does not satisfy the given 
condition. 



6 Two and four modes 

Now we will consider two simple systems, which can be treated explicitly in 
the general case (without using the conditions of theorem 15. ip . The first is 
one Fermionic mode for each party, which is only meant to illustrate that the 
fidelity optimized over quasifree maximally entangled states coincides (at least 
in this simple example) with that over any maximally entangled state (i.e., 
the maximally entangled fraction). The second is the first non-trivial case, two 
Fermionic modes for each party. Here we are not looking at the optimization of 
the fidelity along the lines of Theorem l5.1[ but we postpone, as already proposed 
at the end of Section HJ the twirling in Step 2 and do it only at the end after 
the parity measurement. This leads to a two qubit system, whose maximally 
entangled fraction can be explicitly calculated and therefore compared to the 
fidelity our protocol provides. 

In an appropriately chosen real basis e^ 2 ,e^ 2 (cf. equation (|90p ). the off- 
diagonal blocks of a general covariance matrix S are diagonal (cf. the corre- 
sponding discussion in Appendix IC|) and S therefore becomes 

1 ia ic \ 

ia 1 id , , 

ic 1 ib ( ° 8j 

-id -ib 1 / 

with a, b,c,d e K, and the constraint < S < 1 now reads 

1 + (ab - cdf >a 2 + b 2 + c 2 + d 2 <2. (59) 
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To get the corresponding density matrix it is convenient to use the Jordan- 
Wigner transformation 



V2 



a*®!, B(e B ) 



71' 



i = 1,2. 



(60) 



/ 1 








b \ 








d 








— c 








V « 








ab — cd J 



By means of the Wick theorem we can then determine the correlation matrix 
r ij = Ps(& 1 ® f" 7 ) where er° = I. The matrix turns out to have the form 



(61) 



Now consider a quasifree maximally entangled state tpE- The corresponding 
basis projection is parametrized by the real orthogonal matrix given by equation 
(f2"T)l . which is now 2x2. The fidelity (ipE, Ps^e) maximized over all such E 
then becomes 

sup^B, psiPe) = 7 max{l — ab + cd+ \c + d\, 1 + ab — cd+ \c — d\} (62) 
Y 4 

where the supremum is attained by Y = sign(c + d)I if the first and by Y — 
sign(c — d) diag(— 1, 1) if the second expression is larger. Note that the right 
hand side of (|62j) coincide with the supremum of (ip,pgip) over all (i.e. not 
necessarily quasifree) maximally entangled states [B]. This indicates that we do 
not give away distillation quality by restricting in Step 2 to maximally entangled 
quasifree states. 

If we look at it as a Fermionic bipartite system, the two mode case is too 
simple, because the local algebras A, B (cf. Equations (|20[) and (|21[ll are Abelian 
- hence no quantum degrees of freedom remain locally. The easiest non-trivial 
case requires two mode on each side, hence four modes in total. Unfortunately 
the corresponding covariance matrix S has 16 free parameter (even after diago- 
nalizing the off-diagonal blocks) and it seems to be reasonable to make further 
simplifying assumption. One possible choice is to assume (considering again an 
appropriately chosen real basis) that the submatrix Y$ is not just diagonal, but 
a multiple of the identity, i.e. Ys = crl and that the commutator [Xs, Zg] van- 
ishes. We can then transform the submatrices Xg and Zg to a normal 
and S becomes 



b ~ 2 



1 iv\ 

—%V\ 1 




ia 

ia 


\ 




1 iv 2 
-iv 2 1 




ia 

ia 


— ia 

—io 




1 iv% 
-iv z 1 






—ia 

~ia 




1 1V4 
—iv£ 1 J 



(63) 



7 An antisymmetric matrix can be always be transformed to a form of having only two- 
dimensional antisymmetric blocks in the diagonal by a unitary. When Xg and Zg commute 
this can be done simultaneously by a Bogolubov unitary, obviously leaving Yg = crl invariant. 
In this case, the state is a product of 2-qubit states. 
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In this case basically everything can be calculated explicitly, however, many of 
the results are too complicated and not very useful. Therefore we will only give 
a summary of the most important points. 

• The probability to get the outcome H — h or during a joint parity mea- 
surement is 

l + (-4) d Pf(S-l/2) 1 + (a 2 - ^ 3 )(a 2 - v 2 v A ) 
P = j = g (64) 

• After such a measurement the (untwirled) posteriori states p^~ and pg + 
are always separable, while p^ + and p^~ can be entangled. Their con- 
currence, when non-vanishing, has the form ~ (4er 2 — ^2, ^3, ^4)1) > 
which shows the expected behaviour: a ^ is necessary for entanglement. 

• The entanglement fidelity / = f(S, 2, 1, Vs) from Equation (l4"6"|) where Vs 
is the partial isometry given by the polar decomposition Yg = Vs\Ys\ has 
the form (i.e. this is the output fidelity of the optimal protocol discussed 
in Section [5]) 

(1^3 - o- 2 -l) (^4-ct 2 -1)+4ct 2 
1 = 8^ ■ (65) 

The interesting question is now, whether the "canonical" choice V — Vs for 
the partial isometry V (and therefore the corresponding quasifree, maximally 
entangled state maximizes the quantity pf or what is equivalent in this casq3 
the fidelity /. Direct optimization of the latter (e.g. if the isoclinic decomposition 
of SO (4) is used to parameterize the set of quasifree, maximally entangled states) 
can be done explicitly, but is quite difficult (and therefore omitted here). 

Fortunately, we can easily trace the problem back to a known result. To this 
end recall the remark from the end of Section @] that we can basically swap the 
steps 2 and 3 of our protocol and do the joint parity measurement directly on 
the state ps- If the measured parities coincides, and if we identify the Hilbert 

f2) (2) f2) (2) a 

spaces n ( +> ®HX and U k _>®Kl' as described in SectionH we get the mixture 

P = ~Ps + + ^~p~ p s ~ ' ( 66 ) 

as the output state (cf. Equation (|47p). Twirling with the symmetry group of 
4>e,+ — ipB,— leads to the same isotropic state we would get if we run the 
protocol in the usual order (cf. again Section!?]). This implies in particular that 
we can calculate the fidelity (ipE,±, pipE,±) in terms of Equation (|42j) . i.e. 

KWe,±,PWe,±) = — ■ (67) 

P 

The left hand side of this equation if bounded by the maximal singlet fraction, 
i.e. the supremum of ^) over all maximally entangled two qubit states. 

The latter can be calculated according to [6J which leads to 

supp* m) = max{/, g} , (68) 



8 We do not consider a subsystem, i.e. m = 2 and the unit operator is the only possible 
choice for D. The probability p, however, does not depend on V, but only on D. 
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where / is given in Equation (|65p and g is 

(flU 3 - <7 2 + 1) (p 2 U 4 - CT 2 + 1) 

.9 = ^ • (69) 

Hence we have to determine the parameters for which / is bigger than g, because 
in this case the entanglement of the isotropic output states a ±:iz is as big as it 
could be, and the optimal protocol proposed in Section[S]is really optimal (within 
the given class) although the assumption of Theorem 15.11 (X$ = or Zg = 0) 
is not met. 




-1,4 



Figure 1: The figure shows two surfaces given by the fidelity / of p 4 (x, y, a = 0.2) 
and its maximal singlet fraction max{/, g}. Since / > g for most of the values 
of the parameters, it is only at the left and right side, where we see that the two 
surfaces differ. This region shrinks if we increase a and disappears. 

We have done the last step numerically and chosen x = v% = v\ and 
y = ^4 = v%. In this way we have a product of two identical two- mode state 
p^(x,y,a : a) at hand as seen from Using the notation p 4 (x,y,a) for the 
state, we compare the fidelity / with the maximal entangled fraction (|55)) . The 
result is that / > g for most of the parameter values, and in particular for higher 
fidelities. 

7 Distillation from lattice Fermions 

To get another example for our scheme let us have a look at Fermions propagat- 
ing on an infinite, regular lattice. For simplicity, we will restrict our discussion 
to one dimension, but everything we will say in this section, can be easily gen- 
eralized to higher dimensions. 

The system we are looking at has infinite degrees of freedom. Therefore 
the description used up to now has to be modified a bit (cf [3] [2] for a detailed 
exposition). As a starting point consider the Hilbert space K,^°°^ and the complex 
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conjugation J given by 

Now we can consider as before Majorana operators B(x) smeared by elements 
x G operating irreducibly on a (now infinite dimensional, but separable) 

Hilbert space and satisfying the CAR. As in finite dimensions, a nor- 

malized element Q$ G "H' 00 ) describes a pure, quasifree state if its correlation 
functions are given in terms of Wick's Theorem (cf. Equation ©) and a basis 
projection S G S(/C'°°- ) ), which is defined as in the finite dimensional case. The 
Hilbert space "H^ 00 ' and the operators B(x) can be constructed as in Appendix 
[X] in terms of Fock spaces and creation/annihilation operators, but this is not 
important for us, because the triple consisting of KS°°\ J and S fixes 7^°°) and 
the B(x) uniquely up to unitary equivalence. 

Consider now a finite set A of integers with cardinality n and define the 
projection I A G B(/C(°°)) 

(I A a;)j = xa(j> 3 , jeZ (71) 

where xa denotes the characteristic function of A. This projection commutes 
with J (i.e. it is a real projection) and defines therefore a (finite dimensional) 
subsystem of the infinitely extended system. It is completely described by the 
algebra 

A A = span{S(a;i) • • • B(x t ) \ I G N, I AXj = xflj = l,...,l} (72) 

which is isomorphic to the CAR algebra of n modes. All operations, measure- 
ments, etc., which only affect the modes located in A (and leave the rest of 
the system untouched) can be described in terms of this algebra. It is acting 
reducibly on 'H^ 00 \ but we can cure this defect with a twisted tensor product. 
In analogy to (|17[) we can introduce a unitary 



U A : H {oc) ^ H (n) (g> H {oc) (73) 



such that 



utBvm-l 8 ®** (74) 

w A \0®B(x) xe (I- Ja)/c(°°) s/c(°°), 

where denotes, as before, the Fermionic Fockspace of n modes. Dropping 
the complementary subsystem (which is now infinite dimensional) can be de- 
scribed again by the channel (cf. Equation (|18|) 

B(H in) ) A A (A) = U* A A ® W A g B(H {oo} ). (75) 

After applying this operation, the remaining n modes are in the quasifree state 
ps, a with covariance matrix S A = I A SI A and we end up with the structure 
studied in Section [5J 

Now let us come back to distillation and assume that Alice and Bob share 
an infinite 1-dimensional lattice containing a number (typically infinite) of 
Fermions. They can control only two contiguous regions Aa and A# of L lattice 
sites each and located at a distance of N sites. Everything outside of A = A^UAb 
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can not be manipulated and is therefore ignored^. In other words we apply the 
channel Aa defined in Equation (|75[) and get a bipartite Fermionic system con- 
sisting of 2L modes where the projections I\ A and I\ B play the role of I a and 
1b from Section [31 Hence we can apply the optimal protocol studied in Section 
[5] and calculate the success probability p(L, N) and the fidelity f(L, N) of the 
output systems. Both quantities depend on the geometry of the regions and 
Ab (i-e. their length L and distance N) and to study this functional behavior 
gives a deep insight how distillable entanglement is distributed along the chain. 

Particularly interesting is the question, whether p and / can be made (as a 
function of L) arbitrarily close to 1. This would imply that using a sufficiently 
big chunk of the system is as good as using a sufficiently large number of the 
same system in the same state (which is typically done in ordinary distillation 
protocols such as hashing, and we have assumed this in Section [U too). In this 
sense our result is closely related to the question how much entanglement can 
be distilled with certainty from one copy of an infinitely extended system. This 
is called the systems one copy entanglement and studied in a number of papers 

musing. 



8 Free Fermions hopping on a 1-dimensional lat- 
tice 

As an explicit example for the discussion from the last section let us consider 
the unique ground state of free Fermions hopping on a infinite, 1-dimensional, 
regular latticJ^l (i.e. Z). ft is given by 

where 

1 oo 

3 F^ F(F) e IAS* 1 ) «> C 2 , J 7 (F)(x) = —= V e mx F n (77) 

v 2n . — ' 

j = -oo 

is the Fourier transform and E g B(L 2 (S 1 )) the projection to the upper half- 
circle. Proofs and further references for all of this can be found in [3] . 
As in the last section we will restrict the state to the region A with 

A A = [-L,0), A B = [N,N + L), A = A A UA B . (78) 

This results in the restricted density matrix pi A si A , with a covariance matrix 
IaSIa which becomes in an appropriately chosen real basis 

/ l A /2 + tX iY \ 
^ SI ^{ - lY T 1b/2 + 1 Z ) (79) 

9 We ignore here the fact actions by Alice and Bob on their part of the system usually 
produces disturbances which are propagating along the chain. This approximation is justified 
if the time scale at which Alice and Bob operate is short compared to the propagation speed. 

10 After a Jordan- Wigner transformation this becomes the XX-model. One copy entangle- 
ment of this model is studied in [8]; cf. also the references therein for other related results. 
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with 



and the L x L matrix F r (r = 0, N) is given by 

f if j = k + r 

(Fr)jk = { sm( , (81) 

— A , . , — - otherwise. 

^ [j — k-\-r)7T 

Our goal is now to consider the distillation of qubit pairs from bipartite 
systems in such a state. The crucial step on the algorithmic side is to determine 
the subsystem projection D and the partial isometry V from I\SIa. This is 
done in terms of the singular value decomposition 

2L 

^ = E A ^I^)^I' ^i>^2>...A 2L (82) 

3=1 

of the matrix Y from Equation ([50)1 . This gives rise to 

4 4 4 

Al = I>iWil> D B =J2\^j)(^\, ^ = £|^)<<M- (83) 
3=1 3=1 3=1 

All other quantities (in particular the probability p and the fidelity /) can be 
derived from these three matrices (cf. Theorem 14.11 and 14 . 2[) . 

Unfortunately, even for the matrices -Fo, -FV+l with their very regular struc- 
ture the singular value problem is not explicitly solvable. They are, however, 
Toeplitz matrices and we can at least provide an algorithm to derive D^SD^ 
which is quite efficient (in time and space) and which allows to do numerical 
calculations (at least) in the range L — 10 4 — 10 6 . One possible strategy (and the 
one we have used) is to use iterative procedures like the Lanczos method, which 
only needs matrix-vector multiplications. The latter can then be implemented 
by embedding the Toeplitz matrix T into a bigger circulant matrix C and to 
calculate the product Cx for a column vector x by a discrete Fourier transform. 
The product Ty can then be calculated by first padding y with zeros to fit the 
size of C (which results in a higher dimensional vector x) and projecting out at 
the end the relevant dimensions from Cx. 

Let us discuss now the result of our analysis. Figures [5] and [3] show / as a 
function of L for fixed block distances N = 1, 10, 100, 1000. In Figure [3]we have 
fitted the data with a model function of the form 

f( L ^) = 1 -^m- ( 84 ) 

It is easy to see that for small values of the block length L this fit is not very 
good (in particular for larger block distance N). Therefore we have used only 
data (L, f(N, L)) with L > 20000 for calculating the fit parameters a and b 
(using a least square fit). This difference between the behavior for small and 
large L becomes much clearer in Figure 2] where we have plotted — log(l — /) 
against log(i). The linear part for bigger values of log(i) correspond to the 
1 — b/L a behavior in Figure [31 
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Figure 2: Plot of the entanglement fidelity / as a function of the block length 
L for fixed block distances N = 0, 10, 100, 1000. 




Figure 3: Entanglement fidelity f{L) for N = 1, 10, 100, 1000 fitted with a model 
function (black curves) of the type 1 — b/L a . 
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Figure 4: — log(l — /) as a function of log(L) for N = 1000. For large L the 
functional behavior seems to be linear (cf. the black curve). 
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Figure 5: Fit parameter a as a function of N. The black curve represents the 
average of the values for N > 400. 

21 





i i i i i i i i i i i i i i i i i i i i i N 
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Figure 6: Fit parameter b as a function N. For N > 500 the data is fitted 
logarithmically. 

To get an idea how the two fit parameter a and b depend on N we have made 
a least square fit for the data (L, f(L, N)) with L = 1 + lOOOfc, k g [20, 100) 
and for all N = 0+ lOn with n E [0, 100). The results are shown in Figures[5]and 
[51 We can see that for larger values of N the parameter a is almost constant, 
while b seems to grow sub-linearly. Note that it is dangerous to say more at this 
point, because small fluctuations can be a consequence of the fitting algorithm 
rather than a real physical effect. 

More insight we can get from a look on the function IAN, x) which provides 
the minimal length of the intervals A^/ B we need to geil 1 1 l f(L(N,x),N) > x. 
The result is shown in Figure [7] We can see that the behavior looks linear and 
close inspection of the numerical data indeed shows that at least for large values 
of N the behavior is as linear as it can be for a function which only takes integer 
values. 

Another important function we can investigate is the probability p to get 
during a joint parity measurement on the restricted system in the state pd a sd a 
the results ++ or . 

In Figure [5] we can see that the behavior of p as a function of L looks quite 
similar as f(L) from Figure [3] Again we can find a good fit for larger L with a 
function of the type 1 — a/L b . This is also confirmed by Figure [9] where we have 
again plotted — log(l — p) as function of \og(p). The behavior of p for small L 
and N seems, however, to be more irregular than the one for / and if we look 
very closely at Figure [5] we can see that the linear fit is not as good as in Figure 
U] (the blue curve is bending a little bit around the black line) . 

11 Note that this is the quantity proposed to be studied in I13| . 
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Figure 7: Minimal interval length L(N,x) needed to get a fidelity / > x if the 
distance is N. 




Figure 8: Probability p(L, N) for N = 1, 10, 100, 1000 fitted with a model func- 
tion (black curves) of the type 1 — b/L a . 
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Figure 9: — log(l — p) as a function of log(L) for N — 1000. For large L the 
functional behavior is linearly approximated. 

9 Conclusions 

We have presented a distillation scheme which can be applied to a bipartite 
Fermionic system in a generic quasifree state, and while it is in general not 
optimal we have found substantial evidence that we do not waste too much 
entanglement originally hidden in the given state. The asymptotic rate of this 
protocol (or equivalently the entanglement fidelity of the isotropic states we get 
at an intermediate step) can be calculated explicitly and can therefore be used 
as tool to analyze entanglement of Fermionic systems. 

The latter is particularly true for lattice Fermions. The analysis of the pre- 
vious section showed that a numerical study is possible for effectively very large 
systems up to 10 6 lattice sites. We have only studied a very special quasifree 
state here, but the algorithms used can be easily generalized to other translation 
invariant quasifree states: Their covariance operator is a (2 x 2-matrix valued) 
multiplication operator and the matrix I\SI\ of the restricted state (which can 
be easily calculated by Fourier integrals) is again a (block-)Toeplitz operator. 

The results about free Fermions on a lattice, derived in the last section 
are interesting in their own right and indicate some structure in the fidelity 
f(N, L) and the block length L(N, /). Note that the latter is closely related to a 
quantity proposed in |13) to analyze states with infinite one-copy entanglement. 
The numerical study presented so far does not tell us the whole story of course, 
but it gives us some hints what we might expect from a more detailed, analytic 
study. The latter is of course difficult but a deeper discussion of the asymptotic 
behavior of / and L (which is related to the asymptotic behavior of the singular 
value decomposition of the off-diagonal blocks Y) seems to be a realistic goal. 
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A The self-dual formalism 

The purpose of this appendix to connect the self-dual formalism presented in 
Section [5] to the more familiar description in terms of creation and annihilation 
operators cj , c* , j — 1, n operating on Fermionic Fockspace = J : _(C n ). 
They satisfy the canonical anticommutation relations (CAR) 

{ Cj ,c k } = {c*,c* k } = 0, { Cj ,c* k } = 6 jk l (85) 

and they are irreducible, i.e. for any operator A on we have 

[A, Cj] = [A, c%) = Vj, k = 1, . . . , n A = AJ (86) 

for some A € C. The latter implies in particular that the *-algebra generated by 
the Cj,c* k coincides with Z?^™'). 

Alternatively we can introduce the Majorana operators 

{-4= (c„ + c* ) if a = 1, .... n 

■^{ca-n - c a _ n ) if a = n + l,...,2n. 



The B a are selfadjoint and from (|85| we immediately get 

{B a ,B b } = S aib l- 

Irreducibility of the Cj , c* k implies irrducibility of the B a and vice versa. 



Finally we can smear the B a with a "test function" x £ C 2 ™ 

2n 

B(x) = J2 XaB a , x = (x a ) Q=1 ,..., 2n e K, {n) = C 2n (89) 

a=l 

to get the structure we have started with in Section [5J Note that in this case 
we have = C( 2 "> and J is given by complex conjugation in the canonical 
basis ((Jx) a — x a ). It is easy to see that the B{x) are irreducible, linear in x 
and satisfy the CAR in the form ((4]). 

If we start with the smeared operators we can derive the B a by B a = B(e a ) 
if e a denotes the canonical basis in Ky 1 '. Instead of the e a , however, we can use 
any real basis e a , a = 1, . . . , In of K.^ to get an irreducible family of selfadjoint 
operators B a — B(e a ) on H^ n \ which satisfy the CAR and are therefore 
as good to describe the system as the orginal B a . In other words: The smeared 
operators B(x) provide a coordinate free description of the system, while the 
B a (and the corresponding Cj,cjll) are coordinate dependent. 

This distinction is not purely academic, because the Cj, c* k (and therefore the 
B a from which the Cj , c* k can be derived) are accompanied by a distinguished 
state: The Fock vacuum which is given by the condition cjSl = Vfc = 1, . . . , n. 
The smeared operators B(x), however, does not prefer any state - we have to 
choose a real basis first. 
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To make the last statement more precise let us show how the definition 
of Fock states in Section [5] is related to the ordinary Fock vacuum. To this 
end let us consider a basis projection E £ B(KS n ') and introduce orthonormal 
bases (/i^)j=i,...,n an d (/& )k=i,...,n m the range Rani? C and kernel 

ker E C K.^ of E respectively. Obviously 

f^(/P + /P) ifo = l,...,n 

is a real basis and a short calculation shows that the vacuum state fl e Ti.^ 
given by 

B a =B(e a ), c k = ^={B k + iB k+n ), c* k Q = Vfc = 1, . . . ,n (91) 

satisfies 

(n,B(x 1 )B(x 2 )n) = (Jx u Ex 2 ). (92) 

Validity of the relations in {5} follows from Wick's theorem. Therefore f2 £ H'-™- 1 
is the Fock state with covariance operator E. 

B The parity operator 

In this appendix we will present a proof of Theorem 14.11 and Equation (|13p . 
Unless something else is explicitly stated we will use the notations from Section 
[21 The core result is the following proposition. 

Proposition B.l Let e a G KS n ', a — 1, . . . , In be a real basis of the reference 
space Ky 1 ' and define 

6 = 2 n i n B(e 1 )...B(e 2n ); (93) 
then the following statements hold: 

1. 9 is a self adjoint unitary. 

2. 9 implements the parity automorphism 0, i.e. 9B(x)9 = —B{x) for all 

3. For each unitary R £ B(K,^) commuting with J (i.e. each real orthogonal 
transformation of (Kf- n \J)) we have 

9 = det(R)i n 2 n B(Rei) . . . B{Re 2n ). (94) 

Proof. Item [T] 9 is unitary: Due to B(h)* = B(Jh) and Je a = e a all the 
operators B a are selfadjoint, hence 

9* = 2 n (-i) n B(e 2n )---B( ei ) (95) 

and this implies together with {B(e a ), B(e a )} = 2B(e a ) 2 = 1 that 99* = 9*9 = 
I. 
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To see that 9 is also selfadjoint, note first that B(e a ) and B(eb) anticommute 
iff a 7^ b. Now reverse inductively the ordering of the factors in Equation ([55)1 . 
The first two steps lead to 

9 = -i n 2 n B( ei )B(e 2d ) ■ ■ ■ B(e 2 ) = ~i n 2 n B( ei )B(e 2 )B(e 2n ) ■ ■ ■ B(e 3 ) 

(96) 

Hence we pick up an additional factor — 1 at each second step (i.e. while moving 
a B(e a ) with a odd to the a th position). In this way we get in total a factor 
(— l) n , which shows with Equation (|55|) that 9 is selfadjoint as stated. 

Item [21 With a similar argument we can show that for all a = 1, . . . , 2n we 
have 

9B(e a ) = {-lf^Bie^e = -B(e a )9. (97) 

As a selfadjoint unitary 9 satisfies 9 2 = 1. This implies 9B(e a )9 = —B(e a ), 
which completes the proof of item [21 

Item [21 If R is real orthogonal, the basis e a = Re a is again a real basis and 
the reasoning from above applies. Hence 

9 R = i n 2 n B{R, ei )...B{Re 2n ) (98) 

is a selfadjoint unitary and implements the parity automorphism 0. Only the 
operators ±9 have these properties, and therefore we only have to check the 
sign. To this end consider a quasifree state ps with covariance matrix S. To 
complete the proof we have to show that tr(ps9R,) — det(R)tr(ps9). 

According to the properties of the covariance operator S, we can introduce 
a real, antisymmetric matrix S by 

i§ ab = {e a ,(S--SL/2))e b ). (99) 

With the definition of quasifree states we therefore get 

tr(p s e) = 2 n i n ti(psB( ei ) ■ ■ ■ B(e 2n )) (100) 

n 

= 2'Y* ]Tsign(p) 11(^-1) ^(2,)), (101) 
i=i 

n 

= 2 n i n sign(p) J] iS p{2j _ 1)tP{2j) (102) 
j=i 

= 2"(-l) n Pf(5) (103) 

and similarly 

tr(p s fl ) = 2"(-l)"Pf(iTSi?) (104) 

where we have identified in abuse of notation R with its matrix representation 
in the basis e a . The statement now follows from the properties of Pfaffian (i.e. 
Pf(R*SR) = det(#)Pf(,§)). " □ 

Note that we have shown in addition that the expectation value of 9 in a 
quasifree state ps, is given by 

tr(p s e)=2 n (-l) n Pi(S). (105) 

This is very useful in the proof of Theorem 14.11 The only problem is that we 
have decided about the sign of 9 in an arbitrary way. This can be fixed with the 
following Lemma, which uses the notations introduced in Section [3] 
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Lemma B.2 Consider a bipartite, Fermionic system consisting of 2m modes, 
and a quasifree, maximally entangled state ipE with covariance matrix E G 

1. There is a real basis e a G K,^ 2m \ a — 1, . . . , 4m such that 

ei,...,e 2m G e 2m +i, ■ ■ ■ , e 4m G K, B n) (106) 

and 

Y E e 2m+ j=e2m Vj = l,...,2m (107) 
holds, where Ye = —UaEIb as in Equation {24]) • 
If 9 denotes the parity operator given by Equation {94\ l and this basis, we get 

(i> E ,0i> E ) = (-l) m . (108) 

Proof. Since ip E is maximally entangled by assumption, the off-diagonal oper- 
ator Y is a partial isometry. Hence we can construct the basis e a by choosing 
G2m+i, ■ ■ ■ , 64m as an arbitrary real basis of and defining e±, . . . , ei m by 

(|107ll . In this basis E becomes 

2 \ -il 1 
Inserting this into (|105p we get 

WeMe) =2 2m ¥f{E) (110) 

with 

The Pfaffian of E is easy to calculate and gives (_i)m(2m-i) 2 -2m = (•_ 1 )m 2 -2m 
Inserting this in (JTTDJ) we get the result. □ 

Proof of Theorem\Jl\ Since P + = (1+ 0) /2 holds, the probability p = tr(p s P+) 
can be calculated as p = (l+tr(ps0))/2 where the sign of the parity is chosen by 
(tl>E,6ipE) = 1. Hence the statement follows from Equation (j!05[) and Lemma 
1531 □. 



C The optimality proof 

The purpose of this appendix is to provide the proof of the optimality Theorem 
15.11 The crucial assumption is that the diagonal subblocks of the covariance 
operator S satisfy X = or Z = 0. In both cases the quantity pf we have to 
optimize becomes (this follows directly from Equation (14^1) . and properties of 
the Pfaffian and the determinant): 

pf(S, m, D,V) = ^ (\det(D A YD B + V)\ + \det(D A YD B - V)\) . (112) 

Here we have identified in abuse of notation the spaces DaKS 2 " 1 ^ and DbK,^ 2 " 1 ' 
and interpreted the operators DaYDb ±V as operators on DA& 2m \ such that 
the determinants have a chance to lead to a non-zero value. The first step in 
the proof is to optimize with keeping D fixed. 
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Proposition C.l Consider a bipartite Fermionic system consisting of 2m 
modes and a quasifree state ps with covariance matrix S satisfying X$ = 
or Zg — 0. Assume in addition that the off-diagonal block has maximal rank 
2m, then we have 



2m , * 2m 



pf(S, m, I, V) < pf(S, m, I, V) = [] -i-A + H —± (113) 

fe=i fc=i 

where A& denotes the eigenvalues of \Y\ and V is the partial isometry defined in 
Equation i5ty) . 

Proof. Let us choose a real basis e a in K,( 2m ^ such that Y is diagonalized, 
i.e. (ej, Yek+2m) = fijk^k- The matrix elements Vjk — (ej, Vek+2m) of V in the 



same basis form a orthogonal matrix. Hence the determinants in Equation (j 1 1 2[) 
become 

2m 

X 3 ), (114) 

P 3=1 

where the sum is taken over all permutations of 1, ... , 2m. This can be rewritten 
as a polynomial in the A^: 

det(F + V)= Y[\ 3 det(Vz). (115) 

£C{l,...,2m} J'GS 

Ve denotes the matrix we get from V if we remove the j th row and column 
for all j 6 S. Now assume that both determinants have the same sign. Then 
Equation (|112p leads to 

pf(S,m,l,V)=2 n^ det (^)- (116) 

SC(l,...,2m} j'GS 
S even 

Now note that < Xj < 1 for all j and Vs is a submatrix of an orthogonal 
matrix. Hence |detVs| < 1 with equality iff Vs is orthogonal again. The ex- 
pression in (|116[) is maximal iff det Vs = 1 for all S. But this is implies that all 
submatrix Vs are orthogonal. This is only possible if V is diagonal. A similar 
reasoning holds if both determinants in (|112l) have opposite signs (instead of 
(|116H we get a sum over subsets with an odd cardinality). 

This shows that Y has to be diagonal with eigenvalues ±1. The optimal 
value of pf(S, m, I, V) then becomes 

fid(r+,r_) + fid(r_,r+) (117) 

with 



fid(r +J r_) 



n ( i+a ;) n (-!+**) 

ier+ fcer_ 



II ( 1 + A ^) II (!- A ^ ( 118 ) 

jer + fcer_ 



where T + and T_ are disjoint sets of integers satisfying T + UE = {1, . . . , 2m}. 
Also note that the second equation in (I118[) holds, due to < Xj < 1. The 
proposition is proved if we can show that 

fid({l, . . . , 2m}, 0) + fid(0, {1, . . . , 2m}) > fid(T+, T_) + fid(r_ , T + ) (119) 
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holds. To prove this inequality rewrite fid as 



fid(r + ,r_) = fid(r + ,0) £ (-lji^Ai (120) 
= fid(0,r_) ]T n^- ( 121 ) 

7cr + je 7 



This leads to 



fid(r + ,r_) + fid(r_,r + )= £ ((-i)i%d(r+,0) + fid(0,r + )) JJa,-. (122) 

7 cr_ je-y 

The definition of fid shows that 

fid(r+,0) > fid(0,r+) (123) 

holds. Hence each summand in (|122l) belonging to 7 with odd cardinality is 
negative. In contrast to that we get 

fid({l, . . . , 2m}, 0)+fid(0, {1, . . . , 2m}) = ]T (fid(r+, 0)+(-l)Hfid(0, T+)) J] X j 

7Cr_ j£~/ 

(124) 

where all summands are positive. This completes the proof. □ 

Proof of Theorem \5.1\ We are now ready to complete the proof of Theorem 15. II 
Hence let us consider 2L modes in the quasifree ps, an integer m < L and a 
real projection D — Da © Db € IC^ 2m \ If Vdsd denotes the partial isometry 
given by the polar decomposition of DaYDb, we get with Proposition lC.il 

fc=l fc=l sc{i ) ...,2m}ies 

|S| even 

(125) 

where Xj are the eigenvalues of |Yds£>| given in decreasing order. Obviously this 
quantity is optimized, if the Xj are as big as possible. Hence the theorem is 
proved if we can show that each D satisfies 

Xj>Xj Vj = l,...,2m (126) 

with equality only if D — D. Here Ai, . . . , Xil are the eigenvalues of \Y\. 

To get this bound note first that instead of \Y\ and |Ybsu| we can look at 
the eigenvalues of Y*Y and Y^ SD YrjsD- To get the corresponding bound we 
will use the Courant-Fischer theorem (Theorem 4.2.11 of [TU]) which states that 
the fc th highest eigenvalue Xk of a hermitian matrix A is given by 

X k = sup inf (127) 

p x^O (X,X) 
Fx=0 
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where the supremum is taken over all rank k — 1 projections F. Now we get 

(x,Y*Yx) > (D B x,Y*YD B x) 



\X,X) \x,x) 

^ (x,D B Y*D A YD B x) (x,D B Y*{l-D A )YD B x) 

— / \ 1 / \ ^ y ' 

[X, X) [X, X) 

> (*,DbY*D*YD bX ) 
(x, X) 

where we have used the fact that I — D A is a positive operator. This leads to 
, (x,Y*Yx) , (x,D B Y*D A YD B x) 

SUp ml ; ; > SUp ml ; : . (131) 

p x^Q \X,X) p x^O \X,X) 

Fx=0 Fx=0 

Hence the k th highest eigenvalue of Y*Y is greater or equal to the fc th highest 
eigenvalue of D B Y* D A YD B . Since this bound is attained by D = D A (&D B the 
theorem is proved. □ 



Proof of Corollarv \5.SX We show that increasing m by two decreases pf, i.e. 

pf(S, m, D, V) > pf(S, m + 2,D, V). (132) 
To this end let us introduce first the polynomial 

e n^=n^-n^<^.-,At>) ass) 

SC{l,...,2m}]6S 3=1 3 = 1 

|S| odd 

Then we can express pf(S, m + 2, D, V) as 

pf(S, m+2, D, V) = i(p/(5, m, D, V)+pf(S, to, D, V)\ m+1 X m+ 2+gX m+1 +g\ m+2 ) 

(134) 

The statement follows from Equation (|133[) and the fact that A m+ i < 1 and 

X m +2 < l- n 
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